Spin order in the one-dimensional Kondo and Hund lattices 
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We study numerically the one-dimensional Kondo and Hund lattices consisting of localized spins 
interacting antiferro or ferromagnetically with the itinerant electrons, respectively. Using the Density 
Matrix Renormalization Group we find, for both models and in the small coupling regime, the 
existence of new magnetic phases where the local spins order forming ferromagnetic islands coupled 
antiferromagnetically. Furthermore, by increasing the interaction parameter |J| we find that this 
order evolves toward the ferromagnetic regime through a spiral-like phase with longer characteristic 
wave lengths. These results shed new light on the zero temperature magnetic phase diagram for 
these models. 
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The interplay between charge and spin degrees of free- 
dom in strongly correlated systems has triggered enor- 
mous interest in recent years due to the rich variety 
of phases found in a plethora of compounds. Charge 
and spin superstructures with a doping dependent wave- 
vector were found for example in L^-zSr^NiO^ using 
neutron scattering^ and electron diffraction^. Stripe 
formation together with incommensurate spin fluctua- 
tions in High-Tc superconductors can also be regarded 
as a manifestation of similar phenomena™ as well as the 
charge and spin ordering found in many of the doped 
manganese perovskites. Another large group of com- 
pounds, the heavy fermion materials, present various 
types of ground states including antiferromagnetically or- 
dered states, the normal heavy fermion state as well as 
superconducting and insulating phases. Heavy-fermion 
systems and Kondo insulators are typical examples of 
systems in which the interactions between conduction 
electrons and quantum localized spins are essential^. 
Their physical properties result from an antiferromag- 
netic coupling J between these two types of particles, the 
so-called Kondo lattice model (KLM). The corresponding 
Hamiltonian has the well-known form: 

The first term represents the conduction electron hop- 
ping between nearest-neighbor sites, c| CT ( Ci a ) being stan- 
dard creation (annihilation) operators. In the second 
term the exchange interaction J is antiferromagnetic 
(J < 0), and Bi = ~ J2 a ,a' Aa^aa 1 <Ha> , ("W are Pauli 
matrices). 



It is interesting to note that, in recent years, the same 
model Hamiltonian with ferromagnetic coupling (J > 
0) has been considered to contain the basic physics of 
manganites exhibiting the "colossal" magnetoresistance 
effect 6 7 8 9 . In this case, both localized spins and itin- 
erant electrons originate from manganese d-states. The 
system is assumed to contain essentially Mn 4+ ions with 
three localized t 2g orbitals represented as local spins Si 
and additional itinerant electrons in the e g orbital. Due 
to the strong Hund coupling the spin of the e g electron 
is constrained to be parallel to the local spin on that 
site. Hund's rule together with the hopping term give 
rise to the "double-exchange" (DE) interaction that fa- 
vors ferromagnetic ordering of the local spins 16 . In recent 
literature this model is often referred to as the ferromag- 
netic Kondo lattice (FKLM), however to avoid confusion 
with the Kondo model, we will call it Hund model (HM) . 
The DE mechanism only requires that the system is away 
from half filling and is independent of the sign of J^. 

We will study the Hamiltonian (1) for both, antiferro 
and ferromagnetic couplings, considering S — 1/2 lo- 
calized quantum spins. To this end we will use the 
density-matrix renormalization group (DMRG)i£ with 
open boundary conditions for chains of different sizes. 
We implemented the finite version of the DMRG algo- 
rithm reaching chains sizes of 36 sites (the discarded 
weight being less than 10 -5 ). The different phases 
are characterized through the local spin-spin correlation 
functions and its Fourier transform, the following spin 
structure factor: 

%) = xE e<,M (^^) 

i,3 



2 



where L is the number of sites in the system and n the 
number of conduction electrons per site. 

In Fig.l we present the phase diagram of the one- 
dimensional Kondo and Hund models which we pro- 
pose from our numerical results for several commensurate 
fillings (marked with dashed lines), improving previous 
resultsiiii^iiSiiSiiLiS. As can be seen, except for a scale 
factor, both phase diagrams present great similarities. 
The half-filled case n = 1 is pathological in both models 
whose ground state is very different from the n^l case: 
characterized by a spin gap, however with a different be- 
havior as a function of | J| in the Kondo and Hund cases, 
it has been referred to as a "spin-liquid phase"—. For 
the Hund model it scales to the 5=1 chain with a Hal- 
dane gap. The transition to the FM phase in the Kondo 
model coincides with the previous work mentioned. How- 
ever, for the Hund model we find that the border lies at 
slightly larger values of J as compared with Ref. ^] due 
to the larger system sizes considered here. We also in- 
clude the phase separated regime (PS)— 

At low J | a phase qualified as "paramagnetic" in the 
KLM and "incommensurate" in the HM had been identi- 
fied with exact diagonalization and DMRC 1 ?' 14 ^' 1 ^ 2 ?' 21 , 
this phase is however much less understood than the fer- 
romagnetic phase. In these references the local spin-spin 
correlations are calculated and 5 (q) shows a peak at 
the wave number corresponding to 2kp of the conduc- 
tion electrons. However, no scaling was performed and 
the existence of what we call "island phases" (see below) 
was missed. Recently, the existence of a "spin dimer- 
ized" phase has been reported 2 ^ for the Kondo model at 
quarter-filling (n = 1/2), through the order parameter 

D(i) = (§i ■ The spin structure is of the island 

type ... Till ■■■ similar to the one we have identified pre- 
viously for the DE-Superexchange model— Our present 
calculations reproduce the results for D(i). However, 
the absence of a spin gap2£ and our results on the long 
distance spin-spin correlation functions suggest a non- 
spontaneously broken translation symmetry state com- 
patible with the generalization of the LSM theorem— 
We find that similar structures are present for other com- 
mensurate fillings n = 1/4,1/3,2/3 and 3/4 as well. 
However, we do not find that the "island" phases sta- 
ble at small | J| transform directly into the FM phase- 
Instead, they evolve toward the FM regime through an 
intermediate spiral-like phase with a wave vector q that 
changes from q = 2kp — nn (island) to q = (F) as \ J\/t 
increases. A spiral phase has been proposed analytically 
in a semi-classical model showing a similar evolution to- 
ward the FM phase- 
Let us first present the quarter-filled case n = 1/2 
(Fig. 2). For the Hund model and large J, i.e J/t > 13, 
the FM phase is clearly identified with a peak of 5 (q) 
at q = 0, and the total spin (considering only localized 
spins) is St — N a /2, where N s is the number of localized 
spins in the system. Each electron forms a triplet state 
5 = 1 with the localized spins and all localized spins are 
ferromagnetically ordered. For intermediate values of J 
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FIG. 1: Magnetic phase diagram for the KLM and HM. The 
various phases are described in the text. The circles indicate 
the crossover values. 

(6 < J ^ 10) the FM phase gives rise to a "spiral" phase 
(Fig. 3), characterized by two broad peaks located at 
incommensurate values of q which evolve inwards (the 
momentum of the peak grows from q = to larger values 
and another peak at a symmetric point with respect to 
q = 7r moves towards smaller values). When J decreases 
further (J/t < 6) the spin structure transforms into an 
"island" structure with a more defined peak of 5 (q) at 
q = nir = ir/2 as shown in Fig. 2(f) indicating a four- 
site periodicity. This phase remains stable down to very 
small values of J/t, it has total spin St = and zero 
spin gap. 

Similar results arc obtained in the Kondo case (J < 0). 
The island-type phase is conserved for a small Kondo 
coupling (J/t > —1.3) and the ferromagnetic phase is 
recovered for J/t < —1.6 (Fig. 2 (a), (c) and (e) ). In 
the strong Kondo coupling limit the ferromagnetic state 
differs from the Hund case in the sense that now, the 
conduction electrons form singlet states (5 = 0) with 
the localized spins, so that the system contains N s /2 sin- 
glets and N s /2 unpaired localized spins. The unpaired 
spins are ferromagnetically coupled and the total spin is 
now St — N s (l — n)/2 — N s /4. In the ferromagnetic 
phase, our results for the ground-state energy compare 
well with those obtained with the effective Hamiltonian in 
the strong coupling limit (\J\/t — ► oo) 2 ^ i.e. — J/8 — 2t/ix 
in the Hund case and 3 J/8 — t/-K in the Kondo case. 
Triplet states in the Hund case conserve the hopping t 
while singlet states in the Kondo case acquire a reduced 
hopping t/2. 

The weight of the peak of the spin correlation S(q*) 
for the island (and ferromagnetic) phases increases with 
the chain sizes indicating a quasi-long-range-order in this 
case. Instead, for the spiral-like phase, this weight re- 
mains constant, indicating a short range correlation. A 
similar behavior for these phases is found for the other 
fillings considered. 

In order to visualize the order in real space we show in 
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FIG. 2: Spin structure factor S(q) for the quarter-filled case 
(n = 1/2) for the Kondo (left panel) and Hund (right panel) 
models. From top to bottom: the ferromagnetic, "spiral" and 
island phases, appearing as |J| is decreased. The lower panel 
shows the size dependence of the peak weight S(q*) in island 
(full symbols) and "spiral" phases (empty symbols). 



Fig. 3 the spin-spin correlation functions for the "island" 
and "spiral" phases. For the former case, the correlations 
change sign quite abruptly every two sites indicating that 
the structure is effectively of the island type TTJ4 men- 
tioned above with quasi-long range ordering. The nearest 

neighbor correlation (S(i) ■ S(i + 1)^ changes abruptly 
from site to site, with a slight border effect. Here we 
reproduce the results of Ref. 22. In the latter case, in- 
stead, the "spiral" order is clearly identified and differs 
qualitatively from the island order. The nearest neighbor 

correlation (^S(i) ■ S(i + 1)^ remains ferromagnetic over 

the chain, with a small oscillation of period 2 suggesting 
a dominant spiral-like state with a remnant of the island 
order, which diminishes with J. 

Similar "island" phases are clearly evidenced for n = 
1 /3 at intermediate values of the coupling for both mod- 
els (—0.5 < J/t < 1). Typical examples are for J/t = 1 
and J/t = -0.3 (see Fig. 4(c)- (d)). S (q) shows a 
clear peak at q = nix = 7r/3 and the spin correlation 
in real space presents an "island" structure of three fer- 
romagnetic spins coupled antiferromagnetically between 
islands, basically TtTIII (not shown). The ferromagnetic 
phase as discussed above is clearly recovered for J/t < — 1 
and J/t > 5. In the intermediate region (for values of J 
between these intervals) we find again the "spiral" phase 
(Fig. 4 (a)-(b)). 

A similar behavior occurs for n = 2/3 as shown in 
Fig. 4 (e)-(h). The island structure of the type TTITTI is 
clearly identified for — 1 < J/t < 10 in S(q) with a peak 
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FIG. 3: Real-space local spin-spin correlations for n = 1/2 
and the Kondo (left) and Hund (right) models. The qualita- 
tive difference between the "island" (top panel) and "spiral" 
(lower panel) phases can be seen. 



at q = 27r/3. This island structure has a ferrimagnetic 
character— and we suggest that this configuration cor- 
responds to the new "ferromagnetic" phase reported in 
Ref. |H For 15 < J < 21 we observe a "spiral" phase. 
In the Kondo model we find that the island phase trans- 
forms into the FM regime through another intermediate 
phase, in the region — 3 ^$ J < — 1, with double wave 
length (see below). 
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FIG. 4: Spin structure factors for both, Kondo (left) and 
Hund (right) models showing the "spiral" ((a)- (b) and (e)- 
(f)) and island ((c)- (d) and (g)-(h)) phases for different fill- 
ings. 

To illustrate the crossover between the different phases 
we plot in Fig. [51 the behavior of the wave vector where 
the spin structure factor is maximum, q* , as a function of 
J for several fillings. Here it is clearly seen that for small 
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FIG. 5: Evolution of the momentum where the spin structure 
factor is maximum as a function of the interaction strength J 
covering both models for different fillings. 



\J\ the island phases show up with their characteristic 
wave vector q* — 2kp — mi. This phase is stable for 
a certain region until the "spiral" phase takes over for 
larger \J\. The wave vector of the "spiral" phase decreases 
to zero as | J| increases, leading finally to the FM phase at 
sufficiently large values of the interaction parameter. We 
observe that this behavior is qualitatively similar to the 
"unwinding" of the spiral phase towards the FM order 
obtained by Fazekas and Muller-HartmannAE 5 -. For n = 
2/3 in the Kondo case we find that q* remains fixed at 



7r/3 in the region labeled as "polaronic liquid" in Ref. 
281 In this case, from the real-space spin-spin correlations 
it is difficult to distinguish clearly between a "spiral" and 
an island of type TtTlll but with short range order. 

In conclusion we have studied numerically, using the 
DMRG, the Kondo and Hund models for localized spins 
interacting with itinerant electrons. In addition to the 
ferromagnetic phase at large J|, we find the existence 
of island and "spiral" phases within the "paramagnetic 
regime" in these models. Both phases differ qualitatively 
as seen in different correlation functions. The island 
phase has quasi-long range order and zero spin gap. Fur- 
thermore we show how the ground-state evolves from the 
low \J\/t "island" 2kp-pha.se to the FM regime through 
a spiral-like phase at intermediate couplings. By care- 
fully analyzing the finite-size scaling, we conclude that 
all phases obtained are not due to Friedel oscillations of 
the open boundaries. Based on the results for commensu- 
rate fillings we suggest the phase diagram shown in Fig.l. 
We would like to point out that the difference in the scale 
of J between the Kondo and Hund models has to be a 
consequence of the quantum nature of the localized spins, 
so that a complete understanding of the phase diagram 
requires a full quantum description of the Hamiltonian. 
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